function diff = q_g_p_fun(m,param)

    diff = q_0_p_fun(min(m,param.m_u),param).*q_1_fun(min(m,param.m_u),param) ...
        + q_0_fun(min(m,param.m_u),param).*q_1_p_fun(min(m,param.m_u),param);

end